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We propose a mechanism to change the interaction strengths of a two component condensate. It 
is shown that the application of tt/2 pulses allows us to alter the effective interspecies interaction 
strength as well as the effective interaction strength between particles of the same kind. This 
mechanism provides a simple method to transform spatially stable condensates into unstable once 
and vice versa. It also provides a means to store a squeezed spin state by turning off the interaction 
for the internal states and thus allows to gain control over many body entangled states. 
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I. INTRODUCTION 

The experimental achievement of atomic Bose-Einstein 
condensation (BEC) [[[J |[ [| || has stimulated extensive 
theoretical and experimental studies in this area. One of 
the most impressive examples of applications of BEC is to 
use condensates with internal degrees of freedom to gen- 
erate quantum entanglement [pj, |L @, [l(], , which 
is the essential ingredient for many quantum informa- 
tion protocols It has been shown that the coherent 
collisional interactions in BECs allow to generate sub- 
stantial many-particle entanglement in the spin degrees 
of freedom of a two-component condensate [pj during the 
free evolution of the condensates. The whole time evo- 
lution of the internal degrees of freedom is determined 
by the interaction strengths between the condensed par- 
ticles. Thus it is desirable to control these interaction 
strengths by some external means since this opens the 
possibility to engineer many particle entangled states. 

One possibility is to change the atomic interaction po- 
tential directly by applying an external magnetic field 
which changes the scattering length. If one uses Fesh- 
bach resonances this method allows for con- 

siderable changes in the interaction properties of BECs. 
In this paper we will propose another method to exter- 
nally control the interaction strengths using tt/2 pulses. 
This will not directly change the interatomic interaction 
potential but rather implement an effective Hamiltonian 
with an interaction strength depending on external pa- 
rameters which can be adjusted easily. We will use this 
effective Hamiltonian to study the influence of changing 
the interaction strengths on the external wave function. 
Furthermore, we will be able to control the time evolu- 
tion of the internal degrees of freedom which are solely 
determined by the effective interaction strengths. In par- 
ticular we will show how to turn off the Hamiltonian for 
the internal states, i.e., the internal states will evolve like 
those in an ideal gas. 

Let us briefly explain the basic idea for changing the 
interaction strengths of BECs using tt/2 pulses: If we fix 
the spatial mode function of a two component condensate 
the dynamics of these two components is described by the 
Hamiltoinian |S § 



where J z is the z-component of an angular momentum 
operator J = {J x , J y , J z } given by 

J x = -(alb + tfa), 

J y = ^(tfa-a^b), 

I. = ±{a!a-b% (2) 

where a (b) are bosonic destruction operators for parti- 
cles in internal state 1 (2) with a spatial mode function 
<Pi(2)- The parameter x is determined by the interaction 
properties of the two-component condensate. We apply 
7r/2 Raman-laser or microwave pulses to the condensate 
which rotate the spin around the x or y axis by an angle 
7r/2 depending on the phase of the pulse as we will show 
in Sec. II B. The idea is to use a sequence of 7r/2 pulses 
which rotates H mi to create contributions to the Hamil- 
tonian given by and \Jy f° r a time C St while the 
system evolves with \ Jl for a time St. 

Time averaging leads to an effective Hamiltonian which 
can be written as 



™ ff xJi + Cx(Jy + JZ) _ (i 

^int = 1 , = X- 



1 + 2C 



1 + 2C 



(3) 



xJl 



(1) 



The operator J 2 is a constant of motion and the Hamilto- 
nian therefore is equivalent to H^ t = xJ% up to a (time 
dependent) global phase. Thus, the application of the 
7r/2 pulses effectively leads to a change in the interaction 
parameters from x to x = x(l ~ 0/(1 + %0- For ( = 1 
we find the effective Hamiltonian to vanish. The 
internal states then evolve like those in a non-interacting 
gas. 

The paper is organized as follows. In Sec. |H| we will 
introduce the model. After writing down the Hamilto- 
nian of a two component condensate interacting with a 
classical laser or microwave field we will define a specific 
series of pulses applied to the condensate. Then we will 
calculate the time averaged Hamiltonian and determine 
the dependence of the effective interaction strengths on 
parameters of the external field. In Sec. Ill we will study 
possible applications of this effective Hamiltonian. We 
will investigate the influence of changing the interaction 
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FIG. 1: Schematic plot of the sequence of laser pulses fl/Qo 
against time t in arbitrary units, a) Sequence of tt/2 pulses 
Q(t) applied to the two component condensate between nt c 
and (n + l)t c . The angle a of the tt/2 pulses is shown above 
the corresponding pulse. Graph b) shows the whole series of 
pulses which consists of tt/2 pulse sequences shown in a). 



strength on the spatial wave functions of the condensate 
as well as on the evolution of the internal atomic degrees 
of freedom. Sec. |v| is devoted to the discussion of ap- 
proximations and possible imperfections in our model. 
We conclude in Sec. [v| with a discussion and summary of 
our results. 
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Here, {k,l} € {1 ; 2} and ipk = ?/'fc( x ) is a bosonic field 
operator which annihilates a particle at position x in hy- 
pcrfine state \k). The trapping potential for particles 
in state k is denoted by 14 = Vfc(x) and the mass of 
the atoms is m. The interaction strengths are given 
by u = {uix, 1*22)^1,2} where for simplicity we assume 
Uxi = U22 = u for the interaction between atoms in the 
same internal state, and u\2 ^ u as is the case in Na [ p2| . 
Furthermore we assume the trapping potential for the 
different internal states to be equal, i.e., V\ = V% = V. 
We denote the first excitation energy of V by u and the 
size of the single particle ground state in the potential V 
by a Q . 

Hl describes the interaction of the external field with 
the condensate. We assume that this external field is 
used to apply a sequence of tt/2 pulses to the conden- 
sates. While it is turned on the Rabi frequency is con- 
stant O(i) = Slo . The duration of the tt/2 pulses is thus 
given by td — 7r/2f2o which is assumed to be much shorter 
than the time scale determined by the evolution of the 
condensates due to Hk, i.e., fio ^> u>, and £!o 3> Nu/clq 
PH . Next we investigate the time evolution of the system 
while the external field is turned on. 



II. MODEL 



1. Interaction with the external field 



In this section we start with the Hamiltonian of a two 
component condensate interacting with an external field. 
We investigate the effect of n/2 pulses on the condensate 
and specify a specific series of pulses. We show that this 
specific choice of pulses effectively leads to a change in 
the interaction strengths of the condensate. 



A. Hamiltonian 

We consider a two component BEC consisting of N 
atoms in different atomic hyperfine levels 1 and 2 coupled 
by a time t dependent classical field with Rabi frequency 
fi(t) (internal Josephson effect fl(||l7[ H [OJ)- The clas- 
sical field can be realized by either a Raman-laser or by 
a microwave field applied to the condensate. We assume 
VL to be x independent, i.e. there is negligible momen- 
tum transfer to the condensates due to the interaction 
with the classical field. The Hamiltonian of this system is 
given by H = H BEC + H L , where -Hbec = H 1 (u) + H 2 (u) 
and Hl describe the two component condensate and the 
interaction of external field and condensate, respectively 
pOl |1 El. These terms are given by (with Ti = 1) 



H k (u) = d 3 xipl 




While Hl is turned on, i.e., fl ^ it is the dominant 
part of the Hamiltonian H and we neglect contributions 
of Hk- A 7r/2 pulse is characterized by 



Q(t) = \n(t)\e la , 



\Q(t)\ dt = tt/2, (5) 



with a phase a and it implements the following time evo- 
lution for the bosonic field operators in the Heisenberg 
picture 



1p2 



(6) 



The inverse transformation is given by f/t = U a 1 = 

U a+7r • 



2. Spin operators 

To get an intuitive picture of the effect of tt/2 pulses 
on the condensate we define the spin operator J = 
{J x , J y , J z } by 

J x = - I d 3 x(^ 2 + tpl^l), 
1 /.;:<„ ,,,T t, 

y ~ 2 



J y = o / <rx (V>2>1 - W2), 
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(7) 



In the Heisenberg picture these operators are transformed 
by the tt/2 pulses according to 



UttJzUq — Jy, U- v /2J z U„/2 — Jx 



(8) 



Thus the application of Uq, {U^/2) to the condensate cor- 
responds to a rotation of the spin J around the x(y) axis 
by an angle tt/2, respectively. 



3. Series of tt/2 pulses 

We want to consider a specific sequence of pulses ap- 
plied repeatedly to the condensates. One sequence of 
pulses is shown in Fig. [l]a while Fig. |l]b shows the whole 
series of pulses. As can be seen from Fig. [l]a the conden- 
sate first evolves freely for a time St. Then a tt/2 pulse 
rotates the spin instantaneously around the y axis by an 
angle tt/2 and back after a time £ St. Immediately after- 
wards we rotate the spin by it/ 2 around the x axis and 
then back after time ( St. This sequence of tt/2 pulses 
is repeated as shown in Fig. |l|b. Each sequence takes a 
time t c = (1 + 2()St, neglecting the time td needed to 
apply a pulse. 

There are thus four time scales in this model: (i) td the 
duration of a tt/2 pulse, (ii) t c which is the time needed 
for applying a sequence of tt/2 pulses, (iii) £bec = 1/w 
which determines the time scale of the free evolution of 
the condensates and (iv) tint = a^/uN which is the time 
scale set by the interactions between the particles [ p3[ . 
These four time scales are assumed to satisfy the relation 
td <C t c <C £bec ~ t'mt which can easily be achieved by 
an appropriate choice of external parameters. 



4. Time averaged Hamiltonian 

We now want to study the time evolution of the con- 
densates when the pulses specified in the previous section 
II A 3 are applied to the system. The time evolution op- 
erator U M at time t = Mt c , i.e., after applying M pulse 
sequences is given by 



M 



u 



M 



1=1 



(9) 



where 



U 



U e -ur*aat St Uo 



u -ir/2 e U 7r/2 e 



= e -iH et tta_ 



(10) 



To first order in t c we find for the effective Hamiltonian 




FIG. 2: Interaction strength u (solid curve) for particles in 
the same hyperfine level and U12 (dashed curve) for particles 
in different hyperfine levels as a function of £. For £ = we 
assumed that U12 = U12 = 0.7u — 0.7u. 



where u = {un, U22, U12} with uu = U22 = u and 
_ _ u + (u + uia)C 



(11) 



U12 



l + 2( 
U\2 + 2uQ 
1+2C ' 



(12) 



1 + 2C 



The effective interaction strengths u and u\2 appearing 
in the time averaged Hamiltonian H c r depend on the pa- 
rameter £. Experimentally this parameter can easily be 
changed by adjusting laser or microwave parameters. In 
Fig. 2 we show the dependence of the interaction strength 
on the parameter £. The situation £ = corresponds to 
the case of applying no pulses. We have assumed u > U12, 
i.e., a spatially stable two component condensate. By in- 
creasing C we find u decreases and u±2 increases. At 
C = 1 the two effective interaction strengths are crossing. 
This situation corresponds to the case where the internal 
Hamiltonian H^ t = as discussed in Sec. I. 

Qualitatively we can understand how an effective 
change in the interaction strengths arises by looking at 
the time evolution of a condensate with initially all the 
particles in state |fc). If no laser pulses are applied to the 
condensate the phase accumulated due to the interaction 
during a time £St is proportional to u^St. However, if a 
7r/2 pulse is applied the condensate is put into a super- 
position state of the internal states |1) and |2) according 
to Eq. (H). In this case the phase accumulated during 
a time £St is proportional (it + u\2)C t 5t/2. Then a sec- 
ond 7r/2 pulse brings the condensate back to the internal 
state |fc). Effectively the two 7r/2 pulses thus lead to a 
change of the interaction strength. 
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B. Two mode approximation 



where 



In this section we use the two mode approximation for 
deriving the coupled Gross-Pitaevskii equations (GPE) 
for the condensate. Then we investigate the time evolu- 
tion of the internal atomic degrees of freedom. 



1. External degrees of freedom 

We assume that the condensate can be described using 
only one spatial mode function for each component, i.e., 

ipi(x,t) = a (pi(x,t), and ip 2 (x, t) = b <p 2 (x, t) , 

(13) 

with ipk (x, t) the spatial wave function of condensate k 
and a (b) bosonic annihilation operators for particles in 
condensate 1(2). We put this ansatz into the Hamiltonian 
Eq. (|ll|), assume the state of the condensate to be 



ViVi!iV 2 ! 



-|o>, 



(14) 



with N% particles in condensate 1, N 2 particles in con- 
densate 2 and |0) being the vacuum state. Minimizing 
the expression 



4 t - Hcsl 



(15) 



with respect to the wave functions ipi t2 ( x ) we find the 
coupled GPE equations 

*nsr = f-^ + ^ + «^il^il 2 + «i2^2|^ 2 | 2 Vi' 

at \ Ira J 

= (~?- + V + uN 2 \ l p 2 \ 2 +u 12 N 1 \ ( p 1 \ 2 ) i p 2 , 
at \ 2m J 

(16) 

where ip\p = ipi i2 (x,t) and u and u\ 2 are time depen- 
dent. 



2. Internal degrees of freedom 

We want to simplify the model further by assuming 
N% — N 2 = n being of order 0(y/~N) and set <p\ — p 2 = <p 
in the ansatz for the bosonic field operators. Minimizing 
the terms of order O(N) in expression Eq. ( p"5| ) we find 
the GPE 



.dip ( V 2 Ar u + ui2, > 



(17) 



for the wave function. If ip fulfills the above GPE the 
terms of order 0(y/N) vanish in Eq. (|l^) and the time 
evolution of the internal atomic degrees of freedom is 
given by the Hamiltonian (up to a global phase) 



X = X 



1-C 
1 + 2C 



with X = (u-uu) /d 3 x |^| 4 . (19) 



At ( = 1 the parameter % = 0. There will thus be no 
dynamics of the internal degrees of freedom even for N\ ^ 
N 2 . The internal spin operator J is obtained from J by 
using ansatz Eq. ( |l3| ) and is explicitely given in Eq. (Q). 
A physical interpretation of the evolution of the internal 
atomic states was already given in the introduction Sec. I. 

Putting the two mode ansatz Eq. ( |l5| ) into the time 
evolution operator U M we find 



M 



UnX 



1=1 



U-Tr/ 2 e 



-iffintC St 



U w / 2 e 



— iH int St 



(20) 



m 



eff 



XJzi 



(18) 



where U a is obtained from U a by replacing the bosonic 
field operators V ; i(2) w ith the corresponding annihilation 
operators a (b). 



III. APPLICATIONS 

In this section we study possible applications of chang- 
ing the interaction strength between the condensate dis- 
cussed above. First we show how stable condensates can 
be destabilized and vice versa. Then we investigate how 
a squeezed state of a two component condensate can be 
preserved by turning off the interactions for the internal 
states. 



A. (Un)stabilizing a two component condensate 

The properties of multi-component condensates such 
as spin domain formation have been studied extensively 
both experimentally |2(], ^4|, |g[ [27[] and theoretically 
p9| p8| p9| |30| in the last few years for constant inter- 
action strengths. Here, we study the effect of varying £ 
and thus changing the interaction strengths between the 
condensed particles on the condensate wave functions pk 
and the spatial stability. For u > u\ 2 two initially over- 
lapping condensates should remain spatially stable while 
for u < u\2 they separate. 

We solve numerically Eq. (|H3|) in one spatial dimen- 
sion. The trapping potential is assumed to be har- 
monic V(x) = mu) 2 x 2 /2 with ui the trap frequency and 
ao = \J h/ niLo the ground state size. Figure §a(b) shows 
the condensate wave functions v 3 i(2): respectively. We 
change ((t) as shown in Fig. ^p. Initially C(0) = and 
we assume that u > u\2 1 i.e., the two condensates are 
strongly overlapping. When £(t) > 1 the repulsion be- 
tween atoms in different hyperfine states separates the 
condensates in space since then u\ 2 > u (as can be seen 
from Fig. 2). As soon as £(£) becomes smaller than 1 
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FIG. 3: Condensate wave functions ^ipi as a function of time 
and spatial coordinate Xim are shown in graphs a [b)]. Red 
regions indicate large condensate density, blue regions do not 
contain condensed particles, c) Parameter £ as a function of 
time. The parameters are chosen Ni = 5200, N2 = 4800, 
trapping frequency ui = 40Hz, the interaction strength is as- 



sumed to be 2wi2 



3 • 10 coao. 



again the two condensates become overlapping again. 
Note that the time scale for the separation of the two 
condensates depends also on the imbalance of the con- 
densate particle numbers N\ — N 2 . For N\ = N 2 the two 
condensates do not separate for the parameters chosen in 
Fig. |. 



B. Spin squeezed condensate states 

Next we want to consider the behavior of the internal 
degrees of freedom and show that by choosing £ = 1 and 
thus making the internal Hamiltonian _ff; nt = we can 
store a spin squeezed state. 



separable (i.e. entangled) as has been shown e.g. in ||. 
The parameter £ 2 thus characterizes the atomic entan- 
glement, and states with £ 2 < 1 are often referred to as 
"spin squeezed states" 



2. Preserving a spin squeezed condensate state 



We assume an initial state of the form 



(22) 



created by applying a 7r/2 pulse with a = n/2 to a con- 
densate of particles in internal state |1). The evolution of 
this initial state according to the Hamiltonian H?f t with 
constant x has been studied extensively in ^ and leads to 
one axis squeezing as defined in . Initially £ 2 = 1 and 
is then rapidly reduced. After reaching a minimum value 
the entanglement parameter £ 2 increases again. Our aim 
is to control the interaction parameter \ such that after 
£ 2 has reached its minimum value further evolution of the 
system is suppressed. In Fig. ^ we show a comparison of 
the time evolution with the time evolution operator U-^ t , 
and the effective Hamiltonian Hf^ t . The squeezing pa- 
rameter £ 2 (t) is shown in Fig. |]a and the time dependence 
of ((t) is shown in Fig. |]b. As soon as £ 2 has reached 
its minimum value £ goes rapidly towards one and thus 
prohibits further evolution of £ 2 . The squeezing parame- 
ter £ 2 remains at its minimum value which is close to the 
minimum valued = (3/A0 2 / 3 /2 that can be reached by 
one axis squeezing |3l| . Note that the minimum squeez- 
ing parameter that is reached by the evolution according 
to U^ t is smaller than expected from the Hamiltonian 
H c g as long as t c <C £ B ec is fulfilled (cf. Fig. ||). Wc 
find, however, that this difference is always very small 
and vanishes if we further decrease t c compared to the 
value used in Fig. 0. 



The squeezing parameter 



The entanglement properties of the atoms can be ex- 
pressed in terms of the variances and expectation values 
of the angular momentum operators J. Of particular in- 
terest is the squeezing parameter £ 2 defined by || 



mm 

"12.3 (J n2 



N(AJ n 



\Jn 3 



(21) 



where J n = n • J and the are mutually orthogonal 

unit vectors. If £ 2 < 1 the state of the atoms is non- 



IV. DISCUSSION 



There are two different kinds of approximations in our 
scheme. First, the Hamiltonian H c r we use is time av- 
eraged over the duration of a sequence of pulses t c and 
second, we use a two mode description of the two com- 
ponent condensate for describing the dynamics of the 
internal states. We will discuss these two approxima- 
tions separately since they are independent of each other. 
Also, experimentally it is not possible exactly realize 7r/2 
pulses. Therefore, we will also discuss the influence of 
imperfections in the 7r/2 pulses. 




FIG. 4: a) Squeezing parameter £ 2 as a function of time. 
The dotted line shows the minimum squeezing parameter f 2 ^ 
achievable by one axis squeezing as defined in the text. The 
solid curve shows the squeezing parameter obtained by solv- 
ing the Schrodinger equation using the effective Hamiltonian 
H^ t defined in Eq. (ph. The dashed curve shows the squeez- 
ing parameter obtained from the time evolution operator U^t 
with \tc = 5 TO -3 , b) Parameter £ (solid curve) and resulting 
relative interaction strength \ (dashed curve) as functions of 
time. The numerical calculation was done for N = 50. 



A. Approximations 

1. Time averaging 

While Hl is turned on we neglect the free time evo- 
lution of the condensate due to Hbec completely. Typi- 
cally the time evolution due to the applied n/2 pulses will 
take place on a time scale td of a few ns while the typical 
time scale for the free evolution of the condensate £bec 
is on the order of ms. The neglect of the free evolution of 
the condensates during a pulse will lead to an error of the 
order of ^Abec ~ 10~ 4 and is thus well justified. The 
second step in calculating the effective Hamiltonian is to 
average over one sequence of pulses. This will typically 
lead to an error on the order of £ c Abec- In Fig- [| we 
compare the time evolution according to the time aver- 
aged Hamiltonian with the time evolution given by U^ t 
for t c x — 5 x 10~3 and find a very small deviation be- 
tween the two results. 



2. Two mode approximation 

The form of the effective Hamiltonian H c s is equiv- 
alent to the standard form of the Hamiltonian for two 
component BECs. Therefore we expect the same range 
of validity for H c g as for the original Hamiltonian H. 
This also applies to the two mode approximation intro- 
duced in Sec. II B. 
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FIG. 5: Squeezing parameter £ 2 as a function of time. The 
dotted line shows the minimum squeezing parameter 
achievable by one axis squeezing as defined in the text. The 
solid curve shows the squeezing parameter obtained by solv- 
ing the Schrodinger equation using the effective Hamiltonian 
Hf% defined in Eq. (ph. The dashed curve is obtained by 
calculating the time evolution according to the time evolu- 
tion operator U^ t with x^c = 5 • 10 -3 and a 1% error in 
the intensity of the n/2 pulses. The ensemble average over 
R — 2000 realizations is shown. The other parameters are 
equal to those chosen in Fig. ^. 



B. Imperfect n/2 pulses 

It is experimentally possible to adjust the phase a of 
the 7r/2 pulses very precisely, while it is much harder to 
exactly fulfill the integral condition in Eq. (||). Therefore 
we investigate the influence of a violation of this condi- 
tion on the time evolution of our system. We assume a 
random Gaussian error of 1% in the value of the integral 
in Eq. (||) for each pulse applied to the condensate and 
calculate the resulting time evolution in the two mode ap- 
proximation. Figure |B| shows the result for the squeezing 
parameter averaged over R = 2000 different realizations. 
As can be seen from Fig. || an error in the duration and 
intensity of the pulses does not lead to a qualitatively 
different behavior of the system. For some of the realiza- 
tions we obtain a smaller squeezing parameter £ 2 than 
expected from one axis squeezing. In this case the error 
in the n/2 pulses leads to some two axis squeezing which 
yields a smaller squeezing parameter than pure one axis 
squeezing [31 . 



V. CONCLUSIONS 

In this paper we have introduced a method to change 
the interaction strength of a two component condensate 
by n/2 pulses. We have shown that applying a specific 
series of pulses to the condensate leads to an effective 
time averaged Hamiltonian which is of the form of the 
original two component Hamiltonian with an interaction 
strength depending on parameters of the external field. 

As applications of this scheme we have proposed to use 
this Hamiltonian for turning a stable condensate into an 
unstable one and vice versa. We have also shown that it 
is possible to store a spin squeezed state of a condensate 
for, at least in principle, an arbitrarily long time. 
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Finally, we want to point out that the method to 
change the interaction strengths of BECs discussed in 
this paper can experimentally be realized with current 
technology. It is intended to serve as a tool to gain fur- 
ther insight into the properties of BECs as well as to aim 
in engineering many particle entangled states. 
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